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ABSTRACT 

Some of the most successful statistical methods for obtaining distances of plane- 
tary nebulae (PNe) are based on their apparent sizes and radio emission intensities. 
These methods have the advantage of being "extinction-free" and are especially suited 
to be applied to PNe situated at large distances. A similar method, based on the mid- 
infrared (MIR) emission of PNe, would have the advantage of being applicable to the 
large databases created after the various all-sky or Galactic plane infrared surveys, 
such as IRAS, MSX, ISOGAL, GLIMPSE, etc. In this work we propose a statistical 
method to calculate the distance of PNe based on the apparent nebular radius and 
the MIR flux densities. We show that the specific intensity between 8 and 21 micron 
is proportional to the brightness temperature Tf, at 5 GHz. Using MIR flux densities 
at 8, 12, 15 and 21 microns from the MSX survey, we calibrate the distance scale 
with a statistical method by Stanghellini et al. 2008 (SSV). The database used in the 
calibration consisted of 67 Galactic PNe with MSX counterparts and distances deter- 
mined by SSV. We apply the method to a sample of PNe detected at 8 microns in 
the GLIMPSE infrared survey, and determine the distance of a sample of PNe located 
along the Galactic plane and bulge. 
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1 INTRODUCTION 

The problem of determining the distances of Planetary Neb- 
ulae (PNe) has long been considered a cumbersome issue 
because of the numerous complexities related to the evolu- 
tion of these objects. Classical approaches to this problem 
include trigonometric parallaxes (Harris et al. 1997), the 
expansion method (Terzian 1997), spectroscopic distances 
(Mendez et al. 1988), etc. PNe nuclei are far from being con- 
sidered as "standard candles" since they evolve differently 
according to their masses, a parameter that in most cases is 
not easily determined. When the central star is not visible, 
alternative methods based on properties of the nebula itself 
are required. These methods are often treated as "statisti- 
cal" , and are based on nebular properties such as the ionized 
mass (Shklovsky 1956), the evolution of the nebular optical 
depth with time (Kwok 1985, Marten & Schoenberner 1991), 
the electron density (Barlow 1987), etc. 

The first statistical method proposed in the literature. 



attributed to Shklovsky (1956), is based on the assumption 
that the nebular ionized mass is an invariant. Although this 
hypothesis has some limitations (Maciel & Pottasch 1980), 
the method represents a milestone to the distance problem, 
and was used as a first step to calibrate other similar meth- 
ods that followed (Daub 1982, Cahn, Kaler & Stanghellini 
1992). Methods based on radio quantities have several ad- 
vantages over those derived from visual: the interstellar ex- 
tinction does not play a critical role, even in the cases where 
the PNe are severely affected by interstellar extinction; ra- 
dio fiuxes are available for a large number of objects, etc. 
In addition to that, in the last three decades the increasing 
number of observations carried out at high spatial resolu- 
tion, such as the VLA (Zijlstra et al. 1989), has largely im- 
proved our knowledge of young, compact PNe. On the other 
hand, comparisons among some of these "radio" distance 
scales (Zhang 1995; Phillips 2002, 2004) reveal that they of- 
ten exhibit significant biases, showing that the problem is 
still an open issue. 



* E-mail: ortiz@astro.iag.usp.br 



Mid and far-infrared data constitute an interesting al- 
ternative to approach the distance problem, because PNe are 
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Figure 1. Brightness temperature at 5 GHz as a funetion of the specific intensity at 8 (a), 12 (b), 15 (c) and 21 (d) jiva. The straight Une 
represents the best fit with slope 9 log T(,/9 log = 1. Open circles represent PNe without [C/O] data; filled triangles depict O-rich PNe, 
classified a<;cording to their mid-infraxed spectra; filled squares represent PNe with caj:bona<;eous features in their spectra. PNe evolve 
from the top-right to the bottom-left corner of this figure: the radio emission becomes optically thinner whilst the infrared intensity 
decreases with time. 



generally strong infrared sources and their emission in this 
band is little affected by the interstellar extinction. These 
two characteristics allow the detection of PNe severely ob- 
scured by dust, such as PNe hidden in the galactic plane, 
for example. A relationship between the nebular emission 
at infrared and radio wavelengths was established by Zhang 
& Kwok (1993), who studied the correlation between the 
IRAS specific intensity at 60 jum (/eo) and the brightness 
temperature (Tb) at 5 GHz. These two distance-independent 
parameters can be correlated with the nebular size to com- 
pose a new distance scale. Unfortunately IRAS fluxes suffer 
from poor spatial resolution (~ 1') and the severe crowd- 
ness at low Galactic latitudes, limitations that were par- 



tially overcame by the MSX survey (Egan ct al. 1997, 1999), 
thanks to its better spatial resolution (~ 18"). MSX cov- 
ered the Galactic plane between |6| < 5° in 4 photometric 
bands, centred at 8, 12, 15 and 21 /im. Another milestone 
infrared survey was ISOGAL (Omont et al. 2003), a part 
of the ISO mission designed to obtain mid-infrared images 
of the Galactic bulge and specific regions of the Galactic 
plane with spatial resolution better than 6" mainly at 7 and 
15 /im. Although its sensitivity was almost two orders of 
magnitude deeper than IRAS and MSX, ISOGAL suffered 
from its small sky coverage, limited to only ~ 16 square de- 
grees. Nevertheless, the total number of detections, includ- 
ing point and extended sources, amounts to ~ 10^ objects. 
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Table 1. The MSX bands and the coefficients of equation 4. 



Band 


Ai3o(;im) 


aA 






N 


A 


8.28 


-0.1736 


-0.3899 


+0.7960 


59 


C 


12.13 


-0.2728 


-0.2551 


+0.8527 


15 


D 


14.65 


-0.2491 


-0.3291 


+0.9217 


56 


E 


21.34 


-0.2382 


-0.3242 


+1.0342 


52 



More recently, the Galactic plane and bulge were surveyed 
by GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire, Fazio et al. 2004), a project developed as a 
part of the "Spitzor" missiou. GLIMPSE was restricted to 
|Z| < 65° and carried out with the infrared camera "IRAC", 
which obtained images simultaneously in four mid-infrared 
bands, centred at 3.6, 4.5, 5.8 and 8.0 /im and with spatial 
resolution better than 2". These surveys contain hundreds 
of infrared counterparts of PNe that can be used as an ef- 
ficient tool to derive the distance of PNe severely affected 
by interstellar extinction, especially those situated along the 
Galactic plane. 

The main purpose of the present study is to derive a new 
distance scale based on a minimal set of data: mid-infrared 
intensities and the angular size. Although the present cal- 
ibration is based on MSX flux densities (at 8, 12, 15 and 
21 lira), it can be extended to other existing data extracted 
from past (such as IRAS, ISOGAL, GLIMPSE) or future 
surveys carried out at similar wavelengths. 

This paper is organized as follows: Sect. 2 contains a 
general discussion of the radio and infrared emission of PNe 
and the definition of basic concepts used in the following sec- 
tions; Sect. 3 presents the formulation of the method and its 
calibration; Sect. 4 is devoted to show the results, including 
an application of the method to PNe detected by Spitzer, 
our conclusions are presented in Sect. 5. 



2 THE RADIO AND INFRARED EMISSION 
OF PLANETARY NEBULAE 

Since the formation of a planetary nebula succeeds the 
Asymptotic Giant Branch (AGB), much of the circumstel- 
lar dust formed during this latter phase survives after the 
onset of the ionization of the nebula. The first detection of 
infrared radiation of PN origin is credited to Gillet et al. 
(1967) and since then, hundreds of infrared counterparts of 
PNe have been detected. 

The main characteristics of nebular spectra in the mid- 
infrared arc: (i) the continuum, that results from the thermal 
emission of heated grains in the nebula; it appears stronger 
in newly formed PNe and decreases in intensity as the neb- 
ula evolves, as shown in Fig. 1 (ii) atomic emission lines, 
emitted during transitions between fine structure atomic 
levels; a few of them are strong enough to affect broad- 
band magnitudes between 8 and 25 /im, the most prominent 
are: [Arlll] 9.0 /^m, [SIV] 10.53 ^m, [Ncll] 12.8 urn, [NeV] 
14.3 Atm, [Nelll] 15.4 /um, [Slll] 18.7 /xm, [Arlll] 21.8 /im and 
[NeV] 24.1 /tm; (Hi) molecular bands, in emission or ab- 
sorption, are generally caused by vibration transitions of 
molecules present in ice mantles deposited on the grains; 
they appear as broad emission features such as the amor- 
phous silicate band at 9.7 fjxn (due to Si— O stretching), SiG 



at 11.2 /im, the numerous ice and PAH features between 2 
and 20 fim, etc. 

Among the hypotheses adopted in the present method 
is that the infrared specific intensity of a PN is related to 
the evolutionary status of the nebula and does not depend 
on its distance. This was formerly observed by Kwok (1989) 
and Zhang & Kwok (1990), who studied the IRAS 60 fxia 
emission of PNe. The specific intensity can be evaluated by 
the following equation: 

7a = 1.354 X 10*^, (1) 

where /^(MJy sr~^) is the specific intensity, F\ is the in- 
frared flux density (in Jansky), and 6 is the nebular ra- 
dius, in arcseconds. The specific intensity given by the equa- 
tion above represents only an average value across the PNe, 
which in all cases show variations across the surface. An- 
other hypothesis assumed is the spherical symmetry, even 
though it can be valid also for elliptical nebulae since 9^ 
represents to the product of the nebular major and minor 
semi-axes. 

This work is based on the the mid-infrared flux densities 
provided by the MSX Survey. This choice is justified by the 
following facts: (i) the large number of MSX counterparts of 
PNe and the homogeneousness of those data; (ii) the better 
spatial resohition of the MSX survey, compared to IRAS; 
(Hi) the variety of photometric bands (4) between 8 and 21 
fim. The database adopted in this study is the list of MSX 
counterparts of PNe given by Ortiz et al. (2005). The list 
(table 3 of that paper) contains 216 objects, all of them 
smaller than the spatial resolution of the "SPIRIT" camera 
of 18.3". Therefore, all the nebulae included in this study are 
compact from the MSX point of view, and their photometry 
includes the flux density of the whole nebula. 

The radio emission of a PN is generally expressed in 
terms of brightness temperature, which is proportional to 
the radio intensity at 5 GHz, as follows: 

n = 17.67^^, (2) 

where Tt is given in Kelvin, S'sghz in milijansky (mjy) , and 
the apparent radius 9 in arc seconds. Also here the bright- 
ness temperature represents only an average value across 
the nebula, and the same considerations made for equation 
1 must be observed. 

In this work, all Tt data have been extracted from the 
paper of Siodmiak & Tylenda (2001), which represents a 
compilation of 264 objects observed at 1.4 and 5 GHz, all 
of them using the VLA. Thus, the full set of radio data, 
including Ti, and 9, have been obtained with the same set 
of radiotelescopes, a procedure that helps to assure data 
homogeneity. 

Before we begin to devise the distance scale, we 
shall examine a relationship between these two distance- 
independent parameters: Tj, at 5 GHz and the infrared nebu- 
lar specific intensity 7a. According to Zhang & Kwok (1993), 
who studied IRAS counterparts of PNe observed with the 
VLA, these two quantities are well correlated with each 
other. In Fig. 1, the brightness temperature Tt at 5 GHz 
is plotted against the MSX specific intensities at 8, 12, 15 
and 21 fjxa. The data follow a descending track from the 
top-right to the bottom-left corner, as firstly shown by Volk 
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(1992) and Zhang &: Kwok (1993), who used theoretical evo- 
lutionary models of the central star (Schonberner 1981, 1983; 
Blocker & Schonberner 1990) to perform radiative transfer 
calculations of the parameters of the post-AGB envelope and 
the ncwly-formcd nebula. The brightness temperature Tf, de- 
creases with time because of the decreasing optical depth of 
the radiation at 5 GHz as the nebula expands. Therefore, 
nebulae situated near the top-right corner are younger and 
optically thick at 5 GHz, whereas those near the bottom- 
left corner arc older and their continuum radiation emitted 
at 5 GHz is optically thin. One can see that the maximum 
Tf, value (the optically thick case) in our sample is about 

1.0 X 10*K, which corresponds approximately to the typi- 
cal electron temperature of a PN. The sequence ends near 
Tf, ~ 15 K, whereas in Zhang & Kwok (1993) it goes down to 
Tj, ~ 2 X 10~^ K. This important difference is caused by the 
fact that the PNe selected for the present study, taken from 
Ortiz et al. (2005), arc compact and restricted to objects 
with angular diameter below the MSX spatial resolution of 
18.3". This constraint introduces a bias towards younger 
objects, which generally exhibit larger Ti, values. As for the 
central star, little can be inferred from its position in the di- 
agram: Zhang & Kwok (1993) showed that the evolutionary 
track followed by central stars of different masses is practi- 
cally the same, even though there are significant differences 
among the speed that central stars of difiierent masses cross 
the diagram. The correlation between 7^ and Tj, shown in 
Fig. 1 confirms that the MSX intensity is proportional to 
the brightness temperature at 5 GHz over a wide range of 
optical depths and this is expected to be valid for other sur- 
veys carried out at similar wavelenghts. 

2.1 Mid- infrared spectra of PNe: the role of the 
[C/O] ratio 

The mid-infrared spectra of PNe show large differences from 
object to object and the characteristics of a PN are very dif- 
ficult to predict because of the large number of factors that 
play a role in the composition of the spectrum. For exam- 
ple, in a study of ~ 500 PNe carried out by Stasiriska & 
Szcorba (1999) based on IRAS data, a large range of dust- 
to-gas mass ratios was observed, even though this parameter 
seemed not to vary as the nebula evolves. The same study 
also showed that the contamination of the IRAS broad-band 
fluxes by atomic lines is more evident in the [12] and [25] 
bands, especially among the more evolved nebulae. Other 
papers confirm the difficulty in disentangling the individ- 
ual contributions of stellar, nebular and dust contribution 
to the total nebular spectrum (Volk 1992, Zhang & Kwok 
1991). Most of these studies are based on models intended 
to generate synthetic spectra according to a grid of various 
parameters such as the temperature of the central star, elec- 
tron density, dust-to-gas ratio, chemical abundances, etc. 

The effect of the chemical abundances on the final in- 
frared spectrum is also difficult to evaluate because they 
do not only affect the intensities of atomic lines but also 
change the balance of energy in the nebula. Oxygen lines 
for instance act as a cooler of the nebula that eventually de- 
crease the thermal omission of dust. Understanding the role 
of the abundances in the nebular spectrum is also difficult 
because the abundances of certain elements, such as carbon 
for example, are very difficult to determine because of the 



lack of atomic lines in the visible region of the spectrum. 
In addition to that, a few chemical elements show strong 
atomic lines in the mid-infrared, the most important are S, 
Ne and Ar. 

In a study of 25 PNe seen against the Galactic bulge, 
Casassus et al. (2001a) classified PNe according to the 
[C/O] ratio using molecular features identified in mid- 
infrared spectroscopy between 8—13 /xm. Among the 13 
bulge members with identified dust types, 30% were asso- 
ciated with carbonaceous grains, whilst the rcmaiuiiig 70% 
showed SiO features normally Eissociated with O-rich grains. 
However, the same methodology applied to 54 PNe dis- 
tributed throughout the Galactic disc resulted that 30% 
of all chemically identified PNe inside the Solar circle are 
O-rich PNe, but this percentage is only 14% outside it 
(Casassus et al. 2001b). Therefore, the relative [C/O] neb- 
ular abundance shows a radial gradient similar to what is 
observed in the general distribution of carbon and oxygen- 
rich AGE stars (OH/IR, for example, Blorimiacrt et al. 1993, 
Gugfielmo et al. 1998). 

In this study we attempt to use the [C/O] mid-infraxed 
classification to veriiy if the [C/O] ratio produces significant 
deviations from the general evolution of PNe, with respect 
to their broad-band mid-infrared fluxes. Unfortunately, the 
number of PNe with carbonaceous or silicate features iden- 
tified in their mid-infrared spectra is very small, and our 
sample contains only 20 objects classified according to this 
criterium. Among them, only 4 nebulae have their sizes ac- 
curately determined by the VLA: NGC6537, Ml-12, Hb-5 
(all C-rich), and Hl-40 (O-rich). Their position at the up- 
per part of the evolutionary track seen in the Ix versus Tb 
diagram (Fig. 1) indicates that they are among the youngest 
objects in the sample. Actually, the fact that they all show 
a strong continuum is a direct consequence of their early 
evolutionary status. As the PN nucleus (PNN) evolves, the 
circumstcUar dust is slowly destroyed by the UV photons 
originated from the PNN. Therefore, since the observation 
of molecular features in the mid-infrared requires the pres- 
ence of a strong continuum, it is straightforward to conclude 
that any classification based on these features is biased to- 
wards the less evolved PNe. Because of the small number 
of objects with mid-infrared spectra, little can be concluded 
about the role of the [C/O] ratio in the Ix versus Tb lelar 
tionship. 



3 THE CALIBRATION OF A DISTANCE 

SCALE BASED ON MID-INFRARED DATA 

Statistical distances axe generally based on the hypothesis 
that some nebular properties are constant or they can be 
evaluated from distance-independent data. The Shklovsky 
(1956) method for example, assumes that the ionized mass 
of the PN is invariant from object to object, and as a result, 
the distance can be expressed basically as a function of the 
incoming flux Fx and the nebular radius (Cahn & Kaler 
1971, Milne & AUer 1975, Milne 1982, Daub 1982). 

Apart from secondary effects involving physical param- 
eters of the gas, such as electron temperature, the filling 
factor and the abundance of helium, the distance of a PN 
can be expressed by the following general equation: 

D « F^e", (3) 
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Table 2. Distances of the calibrators of Stanghellini et al. (2008, SSV) that appear also in the list of Ortiz et al. (2005), in alphabetical 
order. 



# 


PN 


MSX 


6 




Fl2 


Fl5 


F21 




D12 


Dl5 


D21 


Dave 


Dssv 




name 


name 


(") 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


22 


Apl-12 


003.3253—04.6590 


6.00 


0.15 






5.75 


4.34 






3.99 


4.17 


4.61 


20 


Hb4 


003.1733+02.9275 


3.63 


0.76 




4.88 


6.58 


3.97 




3.68 


4.55 


4.07 


5.04 


178 


Hb5 


359.3566-00.9801 


1.70 


5.22 


11.43 


29.73 


58.04 


3.82 


3.20 


3.01 


3.46 


3.37 


1.69 


82 


Hc2-ll 


259.1519+00.9404 


32.50 


0.24 




1.27 




2.06 




2.50 




2.28 


0.89 


84 


Hc2-15 


261.6163+03.0039 


11.90 


0.13 








3.39 








3.39 


2.17 


98 


He2-85 


300.5892—01.1090 


5.10 


0.44 




2.02 


5.24 


3.81 




4.10 


4.30 


4.07 


3.56 


102 


He2-99 


309.0023-04.2404 


8.50 


0.23 




2.39 


6.48 


3.49 




3.32 


3.46 


3.43 


3.77 


106 


He2-lll 


315.0311-00.3707 


6.00 


0.20 




1.48 




4.10 




4.20 




4.15 


3.51 


112 


He2-123 


323.9552+02.4585 


2.30 


0.35 




1.05 


3.46 


5.42 




6.28 


6.14 


5.95 


5.75 


114 


He2-125 


324.2651+02.5883 


1.50 


0.12 




0.95 


3.46 


7.76 




7.40 


7.06 


7.41 


10.41 


111 


He2-132 


323.1214—02.5639 


8.90 


0.12 




1.36 


2.86 


3.86 




3.77 


4.15 


3.92 


3.43 


116 


Hc2-141 


325.4953-04.0062 


6.90 


0.13 




1.30 




4.19 




4.14 




4.17 


3.47 


118 


Hc2-142 


327.1953-02.2045 


1.80 


3.96 


10.06 


10.63 


23.61 


3.92 


3.27 


3.82 


4.21 


3.80 


7.40 


119 


Hc2-143 


327.8487-01.6702 


2.60 


0.74 


1.44 


4.82 


13.71 


4.54 


5.06 


4.12 


4.25 


4.49 


5.25 


128 


Hc2-152 


333.4300+01.1870 


5.50 


1.13 




6.64 


11.00 


3.15 




2.97 


3.52 


3.21 


3.04 


124 


Hc2-161 


331.5436-02.7880 


5.00 


0.12 




0.93 




4.83 




5.00 




4.92 


5.31 


126 


He2-164 


332.0796-03.3613 


8.00 


0.12 






3.21 


3.99 






4.18 


4.08 


2.79 


138 


IC4637 


345.4789+00.1411 


9.30 








10.91 








2.97 


2.97 


2.40 


23 


IC4673 


003.5529-02.4423 


8.50 


0.33 




1.96 


4.01 


3.28 




3.49 


3.88 


3.55 


3.22 


71 


K3-7 


028.7768+02.7024 


3.15 


0.20 




1.31 


2.58 


5.29 




5.35 


5.95 


5.53 


6.12 


65 


K3-11 


023.8099—01.7905 


1.50 








5.18 








6.41 


6.41 


10.80 


74 


K3-18 


032.0159-03.0364 


2.00 


2.28 


3.75 


4.64 


14.91 


4.14 


4.16 


4.53 


4.54 


4.34 


9.60 


30 


Ml-25 


004.9383+04.9365 


1.60 


0.20 




1.40 


6.17 


6.89 




6.58 


6.02 


6.50 


6.61 


158 


Ml-27 


356.5305-02.3940 


4.00 


0.33 


2.13 


2.60 


21.75 


4.42 


4.07 


4.17 


3.31 


3.99 


4.61 


175 


Ml-29 


359.1141-01.7195 


3.50 






4.78 


4.75 






3.75 


4.97 


4.36 


4.13 


38 


Ml-31 


006.4554+02.0150 


3.50 






2.86 


8.68 






4.26 


4.31 


4.28 


5.10 


27 


Ml-35 


003.9229-02.3230 


2.20 


0.38 




3.18 


4.82 


5.44 




4.83 


5.76 


5.35 


6.09 


53 


Ml-39 


015.9275+03.3589 


2.00 


0.60 


3.12 


2.21 


22.05 


5.22 


4.38 


5.46 


4.14 


4.80 


6.37 


40 


Ml-41 


006.7679-02.2531 


38.00 


0.50 




3.02 


3.67 


1.71 




1.92 


2.44 


2.02 


0.85 


55 


Ml-46 


016.4510-01.9754 


6.00 


0.26 


1.39 


1.70 


12.57 


3.93 


4.12 


4.06 


3.31 


3.86 


3.58 


52 


Ml-50 


014.6219-04.3797 


2.80 


0.19 




1.19 


3.62 


5.59 




5.70 


5.70 


5.66 


5.98 


61 


Ml-51 


020.9993-01.1252 


7.50 


1.19 


4.61 


6.30 


17.03 


2.77 


2.81 


2.72 


2.87 


2.79 


2.29 


54 


Ml-54 


016.0655-04.3900 


6.50 


0.13 




1.41 




4.29 




4.14 




4.22 


3.81 


63 


Ml-57 


022.1811-02.4085 


4.00 


0.39 


1.51 


3.02 


8.76 


4.29 


4.47 


4.02 


4.12 


4.22 


4.40 


62 


Ml-58 


022.0690-03.1846 


3.20 


0.24 




1.15 




5.09 




5.49 




5.29 


5.32 


67 


Ml-59 


023.9280-02.3440 


2.40 


0.86 


1.59 


4.30 


4.89 


4.56 


5.02 


4.35 


5.58 


4.88 


5.42 


59 


Ml-60 


019.7946-04.5277 


1.25 


0.36 




1.84 


5.74 


6.84 




6.66 


6.64 


6.71 


9.41 


75 
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032.7856-02.0419 


1.35 


0.17 




1.74 


2.95 


7.57 




6.59 


7.58 


7.25 


8.77 


160 
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356.9514+04.5575 


1.35 






0.91 


2.97 






7.74 


7.57 


7.66 


10.75 
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M2-16 


357.4824-03.2479 


1.10 


0.19 




1.25 


3.30 


8.05 




7.65 


7.89 


7.86 


7.04 


163 


M2-18 


357.4094-03.5480 


0.75 
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3.88 


9.88 






8.60 


9.24 


15.74 


6 
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1.50 
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7.18 








7.18 


10.16 


14 
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002.2269—02.7854 


4.40 


0.71 


1.51 


3.13 


8.98 


3.72 


4.37 


3.86 


3.97 


3.98 


4.26 


26 


M2-30 


003.7894—04.6562 


1.75 


0.15 








6.99 








6.99 


9.44 


70 


M2-44 


028.5976+01.6550 


4.00 


1.21 




1.99 




3.52 




4.46 




3.99 


4.76 


69 


M2-45 


027.7030+00.7052 


3.20 


0.69 




4.96 


8.66 


4.24 




3.82 


4.44 


4.17 


4.50 


161 


M3-7 


357.1186+03.6107 


2.35 




- 


1.18 


3.72 




- 


6.05 


6.00 


6.03 


6.32 


166 


M3-10 


358.2412+03.6343 


1.50 


0.18 




1.40 


3.35 


7.20 




6.72 


7.11 


7.01 


9.33 


152 


M3-14 


355.4412-02.4679 


1.40 


0.19 




1.29 


3.39 


7.29 




7.01 


7.25 


7.18 


5.70 


8 


M3-23 


000.9220-04.8521 


6.00 


0.16 








4.25 








4.25 


4.25 


159 


M3-38 


356.9847+04.4410 


0.90 


0.24 




1.47 


4.41 


8.33 




7.86 


7.86 


8.01 


14.06 


2 


M3-43 


000.1231-01.1453 


1.90 


0.20 




1.54 


4.25 


6.41 




6.07 


6.22 


6.24 


8.53 


90 


My60 


283.8976+02.2686 


3.80 


0.18 




1.01 


2.40 


4.99 




5.36 


5.70 


5.35 


4.80 


110 


Mzl 


322.4909-02.6123 


12.90 


0.20 




1.36 




3.05 




3.34 




3.20 


2.30 


121 


Mz2 


329.3843-02.8720 


11.50 


0.16 




1.10 




3.32 




3.65 




3.48 


2.36 


125 


Mz3 


331.7287-01.0105 


12.70 


38.79 


80.40 


93.19 


311.80 


1.23 


1.13 


1.17 


1.21 


1.18 


1.45 


77 


NGC2440 


234.8379+02.4206 


9.00 


1.16 


2.47 


8.75 


16.49 


2.59 


3.18 


2.36 


2.72 


2.71 


1.36 


79 


NGC2452 


243.3792-01.0384 


9.40 


0.22 






3.29 


3.39 






3.94 


3.67 


2.84 


86 


NGC2792 


265.7520+04.1025 


6.50 


0.27 




1.50 


4.77 


3.79 




4.08 


1.06 


3.98 


3.05 


141 


NGC6302 


349.5082+01.0551 


22.30 


13.60 


23.77 


71.78 


261.30 


1.18 


1.36 


1.04 


1.05 


1.16 


0.74 
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Table 2. continued 



# 


PN 


MSX 


e 


Fg 


Fl2 


Fl5 


F21 




D12 


Dl5 


D21 


Dave 


Dssv 




name 


name 


(") 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


43 


NGC6445 


008.0752+03.9075 


16.60 


0.29 




2.95 


4.95 


2.60 




2.53 


2.97 


2.70 


1.38 


24 


NGC6565 


003.5326-04.6223 


4.50 


0.19 




1.16 


3.35 


4.64 




4.90 


4.98 


4.84 


4.66 


49 


NGC6567 


011.7433-00.6498 


4.40 


1.54 


2.65 


4.00 


5.34 


3.26 


3.74 


3.63 


4.49 


3.78 


3.61 


47 


NGC6578 


010.8183-01.8272 


4.30 






7.15 


13.28 






3.17 


3.64 


3.40 


3.64 


87 


PB3 


269.7392-03.6199 


3.50 


0.18 




1.31 


2.81 


5.18 




5.17 


5.63 


5.33 


4.89 


129 


Pcl-6 


336.2217+01.9726 


3.60 






1.07 








5.39 




5.39 


5.38 


19 


Pe2-12 


002.8530-02.2915 


2.50 


0.22 








5.69 








5.69 


12.19 




5 10 15 

D(kpc) by SSV (2008) 




5 10 15 

D(kpc) by SSV (2008) 




2 4 6 8 10 

D(kpc) by SSV (2008) 




5 10 15 

D(kpc) by SSV (2008) 



Figure 2. Distances of PNe determined by the present method plotted against the sample of calibrators compiled by SSV, for 4 the 
MSX bands: 8/im (a), 12 //m (b), 15 //m (c), and 21 jum (d). 
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where F\ is the flux density — or the brightness tempera- 
ture — 9 \s the nebular radius, and the constants a and j3 
depend on the wavelength considered. The empirical lineax 
relationship between the radio and the infrared intensities 
discussed in Sect. 2 allows us to use the MSX flux densities 
to calibrate a new distance scale based on the parameters 
above. Equation (3) can be re-written in the following form: 

log -Dmir = ttA log Fa + bx log + Cx, (4) 

where /Jmir is the distance of the nebula, Fx is the MIR 
flux density, 6 is the angular radius, and the constants a>,, 
h\ and ca depend on the wavelength. 

A comparison among the various statistical distances 
existing in the literature shows that some of them have 
significant biases when compared with each other. These 
effects can be minimized by calibrating the distance scale 
with objects in the Galactic bulge and/or in the Magellanic 
Clouds, which have distances determined more accurately. 
The method proposed by Cahn et al. (1992) was based on 
Shklovsky's (1956) and improved by Daub (1982). In that 
study, 36 Galactic PNe with distances accurately determined 
by several methods was used to calibrate an empirical rela- 
tionship between the distance and the nebular fiux (of H/3 
or at 5 GHz, when available) and 0. Eventually, that scale 
was re-calibrated using the Magellanic Clouds as a distance 
standard (Stanghellini et al. 2008, hereafter SSV) and as a 
result they obtained a method that describes the distance 
as a function of at 5 GHz and d. SSV also produced a 
list of distances of Galactic PNe to be used as calibrators in 
future distance scales. 

The distance scale proposed in the present work was 
calibrated with the catalogue of Galactic PNe distances by 
SSV. Only PNe that appear both in Table 1 of SSV as well as 
in the compilation of MSX counterparts of PNe (Table 3 of 
Ortiz ct al. 2005) have been considered (67 objects, Table 2). 
PNe with bad-quality MSX photometry and/or without size 
measurements have been discarded. In order to obtain the 
coefficients of equation 4, we fitted the SSV distances, MSX 
flux densities and the nebular radii using the least-squared 
method. The nebular radii 6 used in this work have been 
preferentially extracted from Siodmiak & Tylenda (2001), 
obtained with the VLA. In the cases where this measurement 
was not available, 6 was taken from the compilation of SSV. 
Among the 67 PNe used in the calibration, there are only 
26 objects with apparent radii listed both in Siodmiak & 
Tylenda as well as in SSV and, among these, there are only 
4 discrepant determinations of 0: M2-16, M3-14, Hb5, and 
NGC2440. Apart from these few exceptions the agreement 
between the two radii (optical and radio) is excellent. 

Table 1 presents the MSX bands, their isophotal wave- 
length (Egan et al. 1999) and the fitting coefficients obtained 
with the least-squared method, where the distance D must 
be given in kiloparsecs. Fx in Jansky, and 6 in arc seconds. 
Because the flux quality varies from object to object and 
from band to band, the number A'^ of PNe used in each fit is 
not the same. A simple examination of Table 1 reveals that 
the coefficients ax, bx and cx show a non-negligible variation 
with the wavelength. The values relative to 12 /im are the 
most uncertain because of the small number of objects with 
good-quality photometry, caused by the lack of sensitivity 
in this band. If one takes the average value of these coeffi- 



cients, equation 3 would be written with (a) = —0.23 and 
(13) = -0.32. 



4 RESULTS 

In Table 2, we list the original data used in the fitting of 
equation 4 and the various distances, obtained by SSV and 
in this work. The columns are: (1) Ordering number, ac- 
cording to Ortiz et al. (2005); (2) Name; (3) MSX name; 
(4) Angular radius ("); (5-8) MSX flux densities (Jansky); 
(9-12) Distance (in kpc) obtained from the MSX bands cen- 
tred at 8, 12, 15 and 21 /xm, respectively; (13) Average Dis- 
tance (kpc); (14) Distance according to SSV. A comparison 
between Dssv and Dave is shown in Fig. 3d. The points are 
approximately dispersed near the identity line within ±50%. 
The scattering of the points results mainly from two sources 
of errors: (i) intrinsic errors of the present method (including 
errors in the data); (ii) errors in Dssv, estimated as ~30% 
in the original reference. According to the authors, this fig- 
ure is similar to that observed in other distance scales and 
cannot be significantly reduced by using better calibrators. 

Among the factors that contribute to the scattering seen 
in Fig. 3 is the error associated with the nebular radius 6. 
FortunatclJ^ the majority of the 6 measurements used in 
our fit were determined with the VLA, which gives spatial 
resolution superior to any optical measurement obtained by 
ground-based telescopes. The error in the determination of 
F'mir caused by an error in 9 is: 

(¥)„,. HM?. (=) 

where ag is the error bar associated with the measurement 
of the nebular radius. Assuming ae = 1 "- corresponding 
to a typical seeing at the visual wavelengths - the error in 
Dmir calculated with equation 5 and 6a taken from Table 1 
is 30% - 40% for a PN with 9^1". This error drops down 
to a half of this value for a PN with 9 = 2"and so on. 

The errors in the MSX fiux densities are small, between 
3% and 6% for the bands hsted in Table 1 (Egan et al. 1999). 
Thus, this source of error does not play an important role 
in the evaluation of the total error in D, especially when 
compared to the contribution of the error in 9. 

Table 3 contains the distances calculated for PNe 
present in the sample of Ortiz et al. (2005), but absent in 
the list of SSV. The columns are: (1) Ordering number, ac- 
cording to Ortiz et al. (2005); (2) Name; (3) MSX name; 
(4) Angular radius ("); (5-8) MSX fiux densities (Jansky); 
(9-12) Distance (in kpc) obtained from the flux density at 8, 
12, 15 and 21 ^m, respectively; (13) Average Distance (kpc). 
In order to compare our results with others, we merged Ta- 
bles 2 and 3 and plotted the average distance obtained by 
the present method against other values given by methods 
based on radio fluxes, extracted from several references. The 
distance scale proposed by Zhang (1995) was based on a cor- 
relation between Ti, and the nebular radius (an upgrade to 
the Shklovsky method), and calibrated with a sample of 134 
"standard" nebulae. Figure 3 shows that Zhang's distance 
scale is about ~ 20% shorter than ours (and SSV, by exten- 
sion). Phillips (2002) proposed a statistical method similar 
to Zhang's (1995), but calibrated it with a sample of nearby 
{D < 0.7 kpc) PNe. His distance scale was based on a re- 
vised relationship between the surface brightness and the 
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Table 3. Distances of MSX counterparts of PNe listed in Ortiz et al. (2005) but not listed as a calibrator by Stanghellini et al. (2008, 
SSV), in alphabetical order. 



# 


PN 


MSX 


e 


Fs 


Fl2 


Fl5 


F21 


Dg 


D12 


Dl5 


D21 


Dave 




name 


name 


{") 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


1 


xli-lD 


nnn i Qnn i n/i Qvof; 

uuu. if5yy-i-u4.o ( zo 


n nn 

u.yu 


n 1 1^ 
U. lo 




O 1 Q 

z.lo 


R no 

D.uy 


c nQ 
o.yo 




7 1 
/ . IZ 


7 OQ 
/ .ZO 


7 va 


iuC5 


m 1 7 

xli-1 ( 


oOo.ooio-hUo.uyUo 


n /I 

U.40 


n /17 
U.4 / 


1 f^O 

1 .oz 


0. 10 


1 n ofi 
lU.Zo 


Q 7f^ 

y . ( 


7 7Q 

( . ( y 


Sift 
0. 10 


a nc; 
O.UO 


a AA 
0.44 


iOO 


xli-lo 


Qc;7 fiOf^n_i_no finQO 


n 7c; 
U. 1 o 






1 7fi 
1. ( D 


o.yo 






7 OS 

/ .yo 


a c^ft 
0.00 


a 07 

O.Z / 


1 TO 
Li A 


xii-iy 


ooo.yt)e54H-uo.oyyo 


n TCi 
U. / U 


n 1 K 
U.IO 


1 /I c: 
1.40 


n nc: 

u.yo 


A nn 

4.yy 


1 n r\A 
1U.U4 


7 nr; 
/ .UO 


n c;i 

y.ui 


a oa 

O.ZO 


a 70 
0. ( z 


1 71 


TJi on 


oOo.yoOU-hUo.ZZol 


l.OO 


n 1 R 

U.lD 




l.DO 


Q OR 
O.Zu 


7 HQ 
/ .UO 




ft Of^ 

O.ZO 


ft Q/1 

o.y4 


ft 7/1 
0.(4 


1041: 


rii-zo 


oo( .DU/o-|-ui. /yoi 


1 OK 
l.ZO 


n 

U.OO 




1 QC^ 

i.yo 


Q Cil 

0.01 


0.00 




ft 1^7 
0.0 ( 


7 /17 

( .4 r 


ft ft/1 
0.04 


lOl 


rii-oi 


1 Q1 n no qi /17 
ooo.iolU — uz.yi4/ 


n Qc^ 
U.oO 








Q n7 
o.U ( 








1 1 ft/1 

11. 04 


1 1 ft/1 
11. 04 




rll-oO 


QC;c^ 7QOQ nQ /171 Q 
oOO. / oZy — Uo.4/ io 


n c;c 


1.00 


Q Ad 
0.4D 


7 1/1 
/ .14 


^1 77 
ol. / ( 


7 QO 
( .OZ 


Q1 

o.yi 


ft OQ 

O.Zo 


7ft 
0. / 


ft Q1 
O.Ol 




rli-oD 


oOo.Oi 14 — U4.y iy4 


n /in 

U.4U 






1 /I Q1 
14. ol 


OA 7n 
Z4. / U 






K 77 
0. / ( 


K 7a 

0. ( 


ft Oa 

O.Zo 


ioi 


TT1 /in 


Qf^o 7T /in no fioi o 
ooy. i 14U — uz.oyiy 


U.DO 


1 nc 
l.Uo 


l.ol 


A 7/1 
4.(4 


1 a 

10. OO 


7 Qn 
( .oU 


ft 7ft 
0.(0 


ft C^Q 

o.oo 


ft Qf^ 

0.00 


ft 7Q 
0.(0 


1 n 
lU 


xli-4 ( 


nni onc^n nQ nQnfi 

uui.zyou— Uo-Uoyo 


1 oc; 
l.ZO 








4.U0 








1 01 
i .zl 


7 01 
/ .Zl 







nnn 7t i 7_i_n/i 7not^ 

UUU. i 11 / -t-U4. / UZO 


n vc; 
U. 1 


n 1 K 

U. ID 






A "iA 
4.04 


Q fit; 

y.oo 






a Q7 
0.0 / 


m 
y .ui 






UUo. iyoO-l-Uo.411U 


o nn 
z.uu 








Q 1 Q 
O.lo 








ft f^ft 
0.00 


ft f^ft 
0.00 


1 Q 

io 


tro on 
rlz-zU 


nno vnftQ i ni ao/IV 
UUz. i yDo+Ul.Dy4r 


1 vn 
1. / U 


n 1 Q 
U. lo 






Q nQ 
O.Uo 


v on 
( .ZU 






7 nn 
/ .UU 


7 1 n 
/ . lU 


A 1 


rlDD 


nn7 o7c;q i ai qaao 
UU i .Zi OoH-Ul.c544z 


Q nn 
o.UU 


1 no 
i.Uo 




7 QA 

/ .o4 


1 7 ft/1 
1 ( .D4 


/I no 
4.Uz 




Q /I Q 
0.4o 


Q 00 

o.oz 


Q 77 
O.ti 


4 


TT„0 oKn 


nnn 70/in i nQ 0/17q 
UUU. / z4y-|-Uo.z4 1 o 


o c;n 
z.OU 


n 1 o 
U.lz 




1 1 o 
l.lz 




ft Q1 
D.Ol 




ft nn 
O.UU 




ft 1 ft 

0.10 


( o 


xVO-D 


nQi nQc;7_i_n/i i oqq 

Uol.UoO ( -|-U4.1Zoy 


n 

U.oO 


U.OO 


1 Qn 

i.yu 


Q na 
O.Uo 


1011 
IZ.ll 


in AK 

1U.40 


7 ao 

/ .oZ 


a QO 

o.yz 


a QO 
o.oy 


a ao 
o.oy 


78 


A/T1 1 
iVl i- 1 z 


Zoo.oooy — uo.yzuo 


n on 
U.yu 


n on 

U.yu 


o c^n 
Z.OU 


Q no 

O.UZ 


7 OQ 

/ .yy 


fi «A 
0.04 


K 7n 
0. ( U 


ft c;7 
0.0 ( 


ft ao 

O.oZ 


ft /1Q 
0.4O 


1 7Q 


A/TI OK 


Qf^Q 01^07 nn 7onQ 
ooo.yoy / — UU. / zuy 


1 fid 
l.DU 


07 
O.Z / 


10. OO 


OQ oa 
zy.ZD 


1 7n /in 

1 / U.4U 


Q on 
o.yu 


07 
z.y i 


Q no 

o.uy 


7Q 
Z. ( 


Q 1 7 
0.1* 


156 


Ml-30 


355.9023-04.2676 


1.75 


0.16 




1.16 


4.27 


6.92 




6.69 


6.39 


6.66 


17 


Ml-37 


002.6807-03.4694 


1.25 


0.14 






5.83 


8.05 






6.61 


7.33 


44 


Ml-40 


008.3356-01.1042 


2.15 


1.25 


2.21 


8.86 


16.92 


4.47 


4.72 


3.77 


4.30 


4.31 


51 


Ml-45 


012.6066-02.7011 


1.25 


0.13 




1.16 


5.76 


8.13 




7.47 


6.63 


7.41 


142 


M2-4 


349.8008+04.4649 


1.00 


0.19 




1.99 


4.15 


8.37 




7.04 


7.71 


7.70 


25 


M2-14 


003.6194+03.1836 


1.10 


0.29 


1.40 


1.42 


6.48 


7.49 


6.34 


7.42 


6.72 


6.99 


37 


M2-31 


006.0403-03.6209 


2.00 


0.19 




1.73 


3.86 


6.39 




5.80 


6.26 


6.15 


58 


M3-25 


019.7518+03.2739 


0.65 


0.75 


1.78 


5.35 


15.11 


7.77 


6.79 


6.34 


6.52 


6.85 


179 


M3-44 


359.3842-01.8140 


2.00 






1.17 


12.36 






6.39 


4.75 


5.57 


56 


M4-8 


018.9038+03.6892 


0.65 


0.09 






4.41 


11.32 






8.74 


10.03 


57 


M4-10 


019.2348-02.2459 


0.60 


0.11 




1.40 


3.43 


11.22 




9.08 


9.52 


9.94 


46 


NGC6537 


010.0977+00.7396 


2.35 


3.64 


4.50 


18.12 


40.34 


3.58 


3.80 


3.06 


3.40 


3.46 


68 


Pcl-18 


027.3446-02.1124 


0.60 


0.42 




2.20 


8.78 


8.88 




8.12 


7.61 


8.20 


16 


Th3-27 


002.6923+04.2158 


1.60 


0.13 








7.45 








7.45 



nebular radius and, when compared to our results (Fig. 3, 
top-right), his method produces distances shorter by a fac- 
tor ~ 3. A revision of that method was eventually proposed 
by Phillips (2004), that assumes a relationship between the 
nebular luminosity at 5 GHz and the nebular brightness tem- 
perature, T(,. A comparison between our results and the dis- 
tances produced by the latter method (Dpho4) is presented 
in the bottom-left panel of Fig. 3: there seems to be a bias 
associated with Dmir < 6 kpc, but both methods agree 
within ~ 30% beyond 6 kpc. 

4.1 Distances obtained from other mid-infrared 
surveys: Spitzer 

The present method was calibrated with nebular flux den- 
sities of the MSX survey centred at 8.3, 12.1, 14.7 and 21.3 
pm. This calibration (equation 4 and Table 1) can be ex- 
tended to other surveys if their effective wavelength is simi- 
lar to MSX. As an example, we present in this section dis- 
tances of PNe detected by the Spitzer satellite during the 
GLIMPSE Legacy Survey (Fazio ct al. 2004). As mentioned 
in the introduction, the IRAC camera had one band cen- 
tred at A = 8.0 /im, not very different from the MSX band 
'A' centred at A = 8.3 /^m. Differently from "point-source 



surveys", GLIMPSE recorded full images of PNe and the 
flux density integrated over the total nebular solid angle 
can be easily computed and used to produce distances by 
our method. 

Kwok et al. (2008) and Zhang & Kwok (2009) obtained 
integrated flux densities of over 60 PNe with sensitivity su- 
perior to MSX. A comparison between the GLIMPSE and 
MSX fluxes of a sample of 22 PNe showed the existence of 
a bias expressed as: F(IRAC8.0)/F(MSX8.3) = 0.90 ±0.36 
(Zhang & Kwok 2009). We applied this correction factor 
to GLIMPSE flux densities at 8.0 to calculate Dmir 
and listed these results in Table 4. The first set of data, 
from row 1 to 17, refers to PNe situated along the Galac- 
tic disc, between 10° — 60° away from the Galactic centre. 
The second set, from the 18th row on, refers to PNe sit- 
uated within 5° from the Galactic centre, i.e. seen against 
the bulge. A simple inspection of these two sets of objects 
shows that, on average "disc" PNe have smaller distances 
when compared to "bulge" objects. Actually the distances 
of the "disc" PNe show a wide range of values, from 0.6 to 
7.2 kpc. On the other hand, the distances of the "bulge" 
PNe are more homogeneous, from 3.1 to 9.3 kpc. The only 
exception is the object named PNG358.8-00.0 (Terz2022): 
considering its large infrared flux and apparent size, it is 
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Figure 3. Distances of PNe determined by the present metliod plotted against other distance scales by: (a) Zliang (1995), (b) Phillips 
(2002), (c) Phillips (2004), and (d) Stanghellini et al. (2008, SSV). The auxiliary dashed lines drawn on both sides of the identity line 
represent a deviation of 30% and 50% from equality. 



a clear example of a nearby object seen against the bulge. 
Table 4 also shows how deep is the Spitzer survey when com- 
pared to MSX: the faintest flux density at 8 fim in Tables 2 
and 3 is Fs = 0.11 Jansky (object M4-10), whereas Table 4 
contains plenty of nebulae with one-tenth of this value. 



5 CONCLUSIONS 

In this work we proposed a distance scale of planetary neb- 
ulae based on mid-infrared flux densities and the nebular 
apparent size. A study of two distance-independent nebular 
parameters - the brightness temperature at 5 GHz and the 
mid-infrared speciflc intensity I\ - showed that these two 
characteristics are very well correlated. This relationship, as 



well as the similar formulation of the distance scale - D as 
a function of the nebular monochromatic specific intensity 
and its size - allows us to state that the reliability of our 
method is similar to SSV's. In fact, a comparison between 
the nebular distances calculated by both methods showed 
that they produce similar results within 50%. Nevertheless, 
since both methods are statistical, each of them might be ac- 
curate only within the natural limitation imposed by the set 
of hypotheses assumed to be valid for the whole set of PNe, 
but that may vary significantly from object to object. These 
variations include the different effect of the various spec- 
tral features (continuum, atomic lines, molecular bands) on 
the final flux, deviations from sphericity, etc. However, like 
in many statistical methods, particular deviations from the 
average may not introduce significant errors in the distance 



10 R. Ortiz et al. 



Table 4. Distances of PNe observed by GLIMPSE. The original data were extracted from Kwok et al. (2008, first set of data) and 
Zhang & Kwok (2009, second set below). The columns are (1) PNG name; (2) other name, (3) MSX name; (4) integrated flux density 
(Jansky) at 8 /im, corrected (see text); (5) nebular radius, in arc seconds; (6) Dmir (kpc). 
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scale, but only in the error bars. As an example, according 
to SSV, excluding bipolar PNe from the sample would not 
change their distance scale by more than 5%, and a similar 
figure is expected here. On the other hand, variations of this 
kind may affect the distance error bars, and a comparison 
between the distances determined by SSV and our results 
showed that are equal within ±50%. 

Although the effects of the [C/O] ratio - either of the 
gas or on the grain surface on the nebular spectrum are 
well known, we could not detect its effect on the distance 
scale proposed, because of the small number of objects in 
our sample with observed C and/or O molecular features. 



As an example of application of the present method, we 
derived distances of a sample of PNe observed by the Spitzer 
infrared telescope. This distance scale can be extended to 
any other database containing nebular radii and flux den- 
sities near 8, 12, 15, and 21 /im, such as IRAS, ISOGAL, 
and future infrared surveys of this kind. 

ACKNOWLEDGEMENTS 

We thank Professor W.J. Maciel for his collaboration in 
the early phase of this work. This work has been par- 



A distance scale of PNe based on MIR data 11 



tially funded by the Sao Paulo Agency for Science Support 
FAPESP, grant no. 2010/18835-3. This research has made 
use of NASA's Astrophysics Data System and the SIMB AD 
Astronomical Database. 



REFERENCES 

Barlow, M.J., 1987, MNRAS, 227, 161 

Blocker, T. & Schonberner, D., 1990, A&A, 240, Lll 

Blommaert, J.A.D.L., van der Veen, W.E.C.J. & Habing, H.J., 

1993, A&A, 267, 39 
Cahn, J.H. & Kaler, J.B., 1971, ApJS, 22, 319 
Cahn, J.H., Kalcr, J.B., Stanghellini, L., 1992, A&AS, 94, 399 
Casassus, S., Roche, P.F., Aitkcn, D.K. & Smith, C.H., 2001a, 

MNRAS, 327, 744 
Casassus, S., Roche, P.F., Aitken, D.K. & Smith, C.H., 2001b, 

MNRAS, 320, 424 
Daub, C.T., 1982, ApJ, 260, 612 

Egan, M.P., Price, S.D., Shipman, R.F. & Tedesco, E., 1997, 
A&AS, 29, 1294 

Egan, M.P., Price, S.D., Moshir, M.M. et aL, 1999, The Midcourse 
Space Experiment Point Source Catalog Version 1.2 Explana- 
tory Guide, AFGL Tech. Rep. AFRL-VS-TR-1999-1522 

Fazio, G.G., Hora, J.L., Allen, L.E., Ashby, M.L.N., Barmby, 
P., Deutsch, L.K., Huang, J.-S., Kleiner, S., Marengo, M., 
Megeath, S.T. et ah, 2004, ApJS, 154, 10 

Gillett, E.G., Low, F.J. & Stein, W.A., 1967, ApJ, 149, L97 

Guglielmo, F., Le Bertre, T. & Epchtein, N., 1998, A&A, 334, 
609 

Harris, H.C., Dalm, C.G., Monet, D.G. & Pier, J.R., 1997, in lAU 
Symp. 180; Planetary Nebulae, eds. H.J. Habing & H. Lamers 

(Kluwer: Dordrecht), p. 40 
Kwok, S., 1985, ApJ, 290, 568 

Kwok, S., 1989, in lAU Symp. 131, Planetary Nebulae, cd. S. 

Torres- Peimbert (Kluwer: Reidel), 401 
Kwok, S., Zhang, Y., Koning, N., Huang, H.-H. & Churchwell, 

E., 2008, ApJS, 174, 426 
Maciel, W.J. & Pottasch, S.R., 1980, A&A, 88, 1 
Marten, H. & Schoenberner, D., 1991, A&A, 248, 590 
Mcndcz, R.H., Kudritzki, R.P., Hcrrero, A., Husfeld, D. & Groth, 

H.G., 1988, A&A, 190, 113 
Milne, D.K., 1982, MNRAS, 200, 51P 
Milne, D.K. & AUer, L.H., 1975, A&A, 38, 183 
Omont, A., Gilmore, G.F., Alard, C, Aracil, B., August, T., 

Baliyan, K., Beaulieu, S., Begon, S., Bertou, X., Blommaert, 

J.A.D.L. et al., 2003, A&A, 403, 975 
Ortiz, R., Lorenz-Martins, S., Maciel, W.J., Rangel, E.M., 2005, 

A&A, 431, 565 
Phillips, J.P., 2002, ApJS, 139, 199 
Phillips, J.P., 2004, MNRAS, 353, 589 
Schonberner, D., 1981, A&A, 103, 119 
Schonberner, D., 1983, ApJ, 272, 708 
Shklovsky I., 1956, Astr. Zh., 33, 222 
Siodmiak, N. & Tylenda, R., 2001, A&A, 373, 1032 
Stanghellini, L., Shaw, R.A. & Villaver, E., 2008, ApJ, 689, 194 

(SSV) 

Stasinska, G. & Szczerba, R., 1999, A&A, 352, 297 

Terzian, Y., 1997, in lAU Symp. 180: Planetary Nebulae, eds. 

H.J. Habing & H. Lamers (Kluwer: Dordrecht), p. 29 
Volk, K., 1992, ApJS 80, 347 

Zijlstra, A. A., Pottasch, S.R. & Bignell, C, 1989, A&AS, 79, 329 

Zhang, C.Y., 1995, ApJS, 98, 659 

Zhang, C.Y. & Kwok, S., 1990, A&A, 237, 479 

Zhang, C.Y. & Kwok, S., 1993, ApJS, 88, 137 

Zhang, C.Y. & Kwok, S., 1991, A&A, 250, 179 

Zhang, C.Y. & Kwok, S., 2009, ApJ, 706, 252 



